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I.  EXECUTIVE  SUMMARY 

This  document  is  the  final  technical  report  for  work  performed  under  Air  Force 
Office  of  Scientific  Research  (AFOSR)  Grant  F49620-93-1-0473.  The  tide  of  the  project 
was  "Materials  for  High-Power  Laser  Applications:  Point  Defects  in  KTP  and  ZnGeP2 
Crystals"  and  the  Principal  Investigator  was  Larry  E.  Halliburton,  Physics  Department, 
West  Virginia  University.  This  was  a  one-year  research  program  to  use  electron 
paramagnetic  resonance  (EPR)  and  electron-nuclear  double  resonance  (ENDOR) 
techniques  to  identify  and  characterize  point  defects  in  KTP,  ZnGeP2,  and  CdGeAs2 
crystals.  During  the  grant’s  active  period,  from  June  1,  1993  to  May  31,  1994,  four 
major  studies  were  completed  and  published.  These  are: 

1.  Effect  of  Crystal  Growth  on  Ti^^  Centers  in  KTiOPO_^.  M.  P.  Scripsick,  G. 

J.  Edwards,  L.  E.  Halliburton,  R.  F.  Belt,  and  G.  M.  Loiacono,  J.  Appl. 
Phys.  26,  773  (1994). 

2.  Identification  of  a  Radiation-Induced  Hole  Center  in  KTiOPO.;.  G.  J. 

Edwards,  M.  P.  Scripsick,  L.  E.  Halliburton,  and  R.  F.  Belt,  Phys.  Rev.  B 
4S,  6884  (1993). 

3.  Electron  Paramagnetic  Resonance  Study  of  a  Native  Acceptor  in  As-Grown 

ZnGeP2.  M.  H.  Rakowsky,  W.  K.  Kuhn,  W.  J.  Lauderdale,  L.  E.  Halli¬ 
burton,  G.  J.  Edwards,  M.  P.  Scripsick,  P.  G.  Schunemann,  T.  M. Poliak, 

M.  C.  Ohmer,  and  F.  K.  Hopkins,  Appl.  Phys.  Lett.  64,  1615  (1994). 

4.  Electron  Paramagnetic  Resonance  Spectra  in  As-Grown  CdGeAs2.  L.  E. 

Halliburton,  G.  J.  Edwards,  P.  G.  Schunemann,  and  T.  M.  Pollock,  J.  Appl. 
Phys. (scheduled  for  Jan.  1,  1995  issue). 

In  the  first  of  these  studies,  a  series  of  Ti^"^  centers  were  formed  in  hydrothermal- 
ly  grown  and  flux-grown  potassium  titanyl  phosphate  (KTi0P04  or  KTP).  These  3d^ 
defects  (S  =  1/2)  were  created  with  60-kV  x-rays  at  77  K,  and  electron  paramagnetic 
resonance  (EPR)  and  electron-nuclear  double  resonance  (ENDOR)  data  were  taken  below 
30  K.  The  ENDOR  spectra  show  that  the  two  Ti^"^  centers  having  the  largest  concentra¬ 
tions  in  hydrothermally  grown  KTP  have  a  neighboring  proton,  presumably  in  the  form 
of  an  adjacent  OH"  ion.  In  contrast,  ENDOR  spectra  show  that  neither  of  the  two  Ti^ 
centers  having  the  largest  concentrations  in  flux-grown  KTP  have  a  neighboring  proton. 
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These  significant  differences  in  the  Icx^  environment  of  the  Ti^  centers  help  explain 
why  KTP  crystals  have  shown  differing  susceptibilities  to  gray  tracking. 

In  the  second  study,  electron  paramagnetic  resonance  (EPR)  and  electron-nuclear 
double  resonance  (ENDOR)  techniques  were  used  to  determine  the  structure  of  a 
radiation-induced  trapped-hole  center  in  potassium  titanyl  phosphate  (KTi0P04  or  KTP). 
A  single  crystal  of  KTP  was  irradiated  with  60-kV  x-rays  at  77  K  to  produce  the  S  =  1/2 
hole  center.  Data  were  then  taken  near  30  K.  Spin-Hamiltonian  parameters  describing 
the  g  matrix  and  three  ^^P  hyperfine  matrices  were  obtained  from  the  angular  dependence 
of  the  EPR  and  ENDOR  spectra,  respectively.  The  principal  values  of  the  g  matrix  are 
2.0008,  2.0250,  and  2.0422  and  the  principal  values  of  the  largest  phosphorus  interaction 
are  -6.809  MHz,  -8.125  MHz,  and  -3.848  MHz.  An  analysis  of  the  spin-Hamiltonian 
matrices  (i.e.,  the  magnitudes  of  their  principal  values  and  the  directions  of  their  princi¬ 
pal  axes)  leads  to  a  model  wherein  the  hole  is  localized  on  a  bridging  oxygen  ion  between 
two  titanium  ions.  The  hole  is  stabilized  at  this  site  by  a  nearest-neighbor  potassium 
vacancy.  This  defect  is  observed  in  both  hydrothermally  grown  and  flux-grown  KTP, 
and  it  is  stable  only  up  to  approximately  160  K. 

In  the  third  study,  electron  paramagnetic  resonance  (EPR)  was  used  to  investigate 
an  acceptor  in  as-grown  single  crystals  of  ZnGeP2.  The  spectra  are  characterized  by 
equally  spaced  triplets  with  1:2:1  intensity  ratios  representing  hyperfine  interactions 
(varying  from  35  to  55  G  in  magnitude)  with  two  equivalent  phosphorus  nuclei.  Their 
angular  dependence  shows  that  there  are  four  crystallographically  equivalent  orientations 
of  the  defect.  The  principal  values  of  the  g  matrix  are  2.002,  2.021,  and  2.074  and  the 
corresponding  principal  axes,  at  one  of  the  four  sites,  are  the  [Oil],  [  1  (X)],  and  [0 1 1] 
directions,  respectively.  Two  possible  models  are  suggested  for  this  acceptor;  either  a 
zinc  vacancy  (V^n)  or  a  zinc  ion  on  a  germanium  site  (Zn^g).  It  also  is  suggested  that 
the  acceptor  responsible  for  the  EPR  signal  is  the  same  acceptor,  namely  ALl,  that  gives 
rise  to  a  dominant  near-infrared  absorption  band. 

In  the  fourth  study,  electron  paramagnetic  resonance  (EPR)  was  used  to  investi¬ 
gate  point  defects  in  a  single  crystal  of  as-grown  CdGeAs2.  Spectra  taken  at  17  K  with 
the  magnetic  field  parallel  to  the  c  axis  show  a  broad  signal  (~165  gauss  wide)  centered 
on  gg  =  2.018  and  a  sharper  signal  (~56  gauss  wide)  centered  on  g^,  =  2.015.  Together, 
these  signals  represent  an  unpaired-spin  concentration  of  approximately  5  x  10^^  cm*^.  A 
single  broad  signal  (~165  gauss  wide)  centered  on  g^  =  1.995  is  observed  at  17  K  when 


the  magnetic  field  is  along  the  a  axis.  It  is  suggested  that  these  dominant  EPR-active 
defects  in  CdGeAs2  consist  of  either  a  cation  vacancy  or  an  antisite  cation  with  the  un¬ 
paired  spin  shared  by  the  four  neighboring  As  ions.  Additional  weak  EPR  lines  appear 
when  the  temperature  is  increased  (from  17  to  40  K)  and  the  microwave  power  is 
decreased  (from  7  to  2.3  mW). 

The  results  obtained  during  these  four  studies  are  described  in  detail  in  this  final 
technical  report  for  AFOSR  Grant  F49620-93- 1-0473. 


4 


n.  EFFECT  OF  CRYSTAL  GROWTH  ON  Tp+  CENTERS  IN  KTP 


Titanium  ions  are  the  primary  electron  traps  in  potassium  titanyl  phosphate 
(KTi0P04),  better  known  as  KTP.  In  this  lattice,  the  titaniums  normally  occur  as 
closed-shell  Ti'*'*'  ions  (ignoring  covalent  effects)  and  they  are  converted  to  Ti^'*'  ions 
upon  trapping  an  electron.  The  Ti^''’  ions  have  a  3d^  electron  configuration,  thus  result¬ 
ing  in  S  =  1/2  centers  that  are  easily  observed  with  electron  paramagnetic  resonance 
(EPR)  and  electron-nuclear  double  resonance  (ENDOR)  techniques.  In  addition  to  their 
magnetic  properties,  Ti^"*"  centers  typically  have  one  or  more  characteristic  optical 
absorption  bands  in  the  visible  and  a  broad  emission  band  in  the  near  infrared. 

Interest  in  the  Ti^"*’  electron  traps  has  arisen  because  of  their  suggested  role  in  the 
optical  damage  process  in  KTP.*’^  Gray  tracks  can  be  observed  in  many  KTP  crystals  if 
they  are  exposed  to  high-power  lasers  while  used  as  a  frequency-doubler.  The  laser 
power  density  needed  to  induce  gray  tracks  varies  from  sample  to  sample,  and  this  sug¬ 
gests  that  the  gray-tracking  phenomenon  is  dependent  on  the  nature  and  quantity  of  point 
defects  present  in  a  particular  sample.  The  Ti^  centers  with  their  broad  optical  absorp¬ 
tion  bands  in  the  visible  are  obvious  candidates  for  the  defects  involved  in  gray  tracks. 

Few  details  are  known  about  the  Ti^'*'  centers  in  KTP.  Roelofs^  observed  EPR 
spectra  from  four  distinct  Ti^"*"  centers  in  flux-grown  material.  These  defects  were 
produced  by  electric  fields  or  by  annealing  in  a  hydrogen  atmosphere.  Andreev  and 
Efimov^  reported  seeing  these  same  four  centers  after  room-temperature  x-ray  irradiation 
of  a  flux-grown  crystal.  Scripsick  et  al.'*  observed  a  quite  different  Ti^"'’  center  in 
hydrothermally  grown  material  after  an  x-ray  irradiation  at  77  K.  These  Ti^'*’  centers, 
plus  others  not  yet  reported  in  the  literature,  have  different  g  matrices  and  different 
superhyperfme  interactions;  however,  no  information  has  been  provided  about  the  differ¬ 
ences  in  the  environments  that  surround  the  Ti^^  centers. 

In  this  section,  we  describe  the  results  of  an  ENDOR  survey  of  Ti^"*"  centers  in 
hydrothermally  grown  and  flux-grown  KTP  crystals.  This  experimental  technique  uses 
hyperfine  interactions  with  neighboring  nuclei  to  probe  the  surroundings  of  a  paramagne¬ 
tic  defect.^  We  find  significant  differences  in  the  local  environment  of  the  Tp"^  centers, 
and  we  suggest  that  this  information  may  help  explain  why  KTP  crystals  have  differing 
susceptibilities  to  gray-tracking. 
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The  hydrothermally  grown  crystals  were  provided  by  Airtron  (a  Division  of  Lit¬ 
ton)  in  Charlotte,  North  Carolina,  and  the  flux-grown  crystals  were  provided  by  Crystal 
Associates  in  Waldwick,  New  Jersey.  Samples  for  the  EPR  and  ENDOR  experiments 
were  cut  from  the  larger  as-grown  boules.  The  faces  of  these  samples  were  perpendicular 
to  the  three  high-symmetry  directions  (i.e.,  the  a,  b,  and  c  axes)  and  their  dimensions 
were  approximately  3x3x4  mm^. 

The  KTP  samples  were  irradiated  at  77  K  with  x-rays  to  produce  the  Ti^"^  centers. 
Irradiations  lasted  up  to  four  hours  with  the  x-ray  tube  operating  at  60  kV  and  30  mA. 
After  irradiation,  the  sample  was  transferred  without  warming  into  an  Oxford  Instruments 
ESR-900  helium  gas  flow  system  for  the  subsquent  EPR  or  ENDOR  measurements.  A 
Bruker  ESP  300  spectrometer  was  used  for  both  the  EPR  and  ENDOR  measurements. 
This  spectrometer  operated  at  9.4  GHz  with  100-kHz  static  field  modulation  in  the  EPR 
mode  and  with  frequency  modulation  of  the  rf  at  12.5  kHz  in  the  ENDOR  mode. 

A.  Defect  Production  Mechanism 

Although  the  Ti^'*'  centers  described  in  this  section  were  produced  by  ionizing 
radiation  (i.e.,  x-rays),  it  is  expected  that  similar  Ti^"*"  centers  will  be  produced  by  an 
intense  laser  beam,  thus  making  our  results  directly  relevant  to  the  gray  tracking  prob¬ 
lem.®  An  advantage  of  using  x-rays  is  the  increase  in  the  total  number  of  centers  that  can 
be  made  available  for  study,  i.e.,  defect  formation  occurs  in  a  larger  volume  of  crystal. 
Also,  producing  and  studying  the  defects  at  77  K  eliminates  the  significant  thermal  decay 
of  these  centers  that  occurs  at  higher  temperatures. 

Our  expectation  that  x-rays  and  intense  laser  beams  create  the  same  types  of  Ti^'*' 
centers  in  KTP  comes  from  the  following  arguments.  Both  methods  of  defect  production 
(x-rays  or  lasers)  begin  with  the  formation  of  electron-hole  pairs  in  the  otherwise  perfect 
lattice.  Many  of  the  electrons  and  holes  will  immediately  recombine,  either  radiatively  or 
non-radiatively,  to  restore  the  perfect  lattice;  however,  a  portion  of  the  electrons  and 
holes  will  migrate  sufficiently  far  from  each  other  and  encounter  a  stabilizing  entity, 
either  a  vacancy  or  an  extrinsic  impurity.  Thus,  a  few  of  the  electrons  and  holes  will 
become  stabilized  at  widely  separated  sites  within  the  crystal  and  form  the  point  defects 
that  contribute  to  an  induced  visible  optical  absorption  (e.g.,  gray  tracking).  In  general, 
the  holes  will  be  localized  on  an  oxygen  ion  with  a  neighboring  entity,  such  as  a  potas¬ 
sium  vacancy,  providing  the  stabilizing  force. ^  The  electrons,  on  the  other  hand,  will  be 
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localized  on  a  titanium  ion  with  a  neighboring  entity,  such  as  an  OH"  ion  or  an  oxygen 
vacancy,  providing  the  stabilizing  force. 

The  electron-hole  pair  needed  to  initiate  damage  will  be  created  in  the  bulk  of  the 
crystal  whenever  above-band-gap  photons  are  present.  Ionizing  photons  are  obviously 
present  in  the  case  of  x-ray  irradiation,  but  there  is  some  uncertainty  as  to  how  a  near- 
infrared  laser,  such  as  a  1.06-|im  beam  from  a  Nd-YAG  laser,  would  produce  photons 
with  energies  greater  than  the  band  gap.  The  most  likely  scenario  involves  nonlinear 
processes  to  achieve  wavelengths  at  or  above  the  band  gap  (for  example,  tripling  the  fun¬ 
damental  beam  or  sum-frequency  generation  with  a  fundamental  and  a  doubled  beam). 
Either  of  these  processes  allow  a  1.064-pm  beam  to  form  355-nm  photons  which,  coinci- 
dently,  nearly  match  KTP's  room-temperature  band  edge  at  350  nm.®  It  is  important  to 
note  that  these  nonlinear  processes,  when  associated  with  electron-hole  production,  do  not 
require  phase  matching  conditions  to  be  satisfied.  An  above-band-gap  photon  would 
immediately  create  an  electron-hole  pair  very  near  the  site  where  the  photon  was  formed 
(or  it  will  be  absorbed  via  other  mechanisms)  and  would  not  propagate  significant 
distances  through  the  lattice  or  join  with  other  photons  to  form  a  coherent  beam. 

B.  EPR  and  ENDOR  Results 

In  this  investigation,  a  series  of  Ti^"*"  centers  were  produced  in  hydrothermally 
grown  KTP  crystals  by  x-irradiations  at  77  K.  Following  the  example  of  Roelofs,  we 
label  these  centers  as  A,  B,  etc.  and  add  a  subscript  hyd  to  identify  the  center  as  appear¬ 
ing  in  hydrothermally  grown  crystals.  Figure  1  shows  representative  EPR  spectra  from 
the  two  centers  having  the  greatest  intensity  after  irradiation.  The  spectrum  in  Fig.  1(a) 
was  taken  at  28  K  with  the  magnetic  field  parallel  to  the  b  axis  of  the  crystal.  It  is  center¬ 
ed  on  3695.24  G  and  has  a  gb  value  of  1 .8321 .  Except  for  a  different  orientation  of 
magnetic  field,  this  is  the  same  center  described  in  the  earlier  EPR  study  by  Scripsick  et 
al.'*  We  now  refer  to  this  center  as  Center  A^^.  The  resolved  hyperfine  pattern  associ¬ 
ated  with  this  center  is  explained  by  interactions  with  four  I  =  1/2  nuclei.  Application  of 
the  ENDOR  technique  provides  precise  values  for  these  hyperfine  splittings  and,  most 
important,  it  unambiguously  identifies  the  responsible  nuclei.  The  ENDOR  spectrum  in 
Fig.  2  was  taken  at  15  K  while  monitoring  the  EPR  spectrum  in  Fig.  1(a).  It  clearly 
identifies  the  responsible  nuclei  in  Center  Af^  as  three  phosphorus  and  one  hydrogen  and 
it  gives  first-order  b-axis  hyperfine  splittings  of  14.68  MHz  for  the  hydrogen  nucleus  and 
11.23  MHz,  17.58  MHz,  and  29.02  MHz  for  the  three  phosphorus  nuclei.  Additional 
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pairs  of  ENDOR  lines  with  spittings  of  2.65  MHz  and  0.61  MHz  represent  phosphorus 
hyperfine  interactions  that  are  not  resolved  in  the  EPR  spectrum. 

A  second  and  less  intense  Ti^"*"  center,  present  in  hydrothermally  grown  KTP 
after  irradiation  at  77  K,  is  shown  in  Fig.  1(b).  We  refer  to  this  center  as  Center  B^. 
These  data  were  taken  at  28  K  with  the  magnetic  field  parallel  to  the  crystal's  a  axis. 

This  center  has  a  g,  value  of  1.8172  and  has  eight  resolved  hyperfine  lines  in  its  a-axis 
EPR  spectrum.  These  lines  represent  interactions  with  three  neighboring  I  =  1/2  nuclei 
which  the  ENDOR  technique  identifies  as  one  hydrogen  and  two  phosphorus.  The  a-axis 
hyperfine  splittings  for  this  center,  obtained  from  the  ENDOR  spectrum,  are  11.00  MHz 
for  the  hydrogen  and  22.88  MHz  and  27.89  MHz  for  the  two  phosphorus.  In  addition  to 
these  three  nuclei,  the  ENDOR  spectrum  shows  that  several  distant  hydrogens  (approxi¬ 
mately  7  A  or  more  away)  are  weakly  interacting  with  the  Ti^  center. 

A  different  series  of  Ti^'*'  centers  were  formed  in  flux-grown  KTP  during  x-ray 
irradiations  at  77  K.  The  Ti^"^  defect  having  the  largest  concentration  corresponds  to 
"Center  B"  in  Roelofs'  paper*  and  we  expand  his  notation  by  adding  a  subscript^  to 
denote  that  this  center  appears  in  flux-grown  crystals.  Thus  we  refer  to  this  center  as 
Center  By^^.  The  EPR  spectrum  from  this  center,  taken  at  25  K  with  the  magnetic  field 
along  the  b  axis,  is  shown  in  Fig.  3(a).  Its  gj,  value  is  1.8614  and  there  are  resolved  hy¬ 
perfine  interactions  with  three  neighboring  I  =  1/2  nuclei,  two  of  which  are  degenerate. 
The  corresponding  ENDOR  spectrum,  taken  at  15  K  and  shown  in  Fig.  4,  identifies  these 
three  nuclei  as  phosphorus.  More  importantly,  there  are  no  ENDOR  lines  which  can  be 
assigned  to  hydrogen.  The  b-axis  hyperfine  splittings,  from  the  ENDOR  spectrum,  are 
9. 10  MHz  for  one  phosphorus  and  23.13  MHz  for  the  other  two  phosphorus  nuclei. 

A  second  radiation-induced  Ti^'*’  center  in  the  flux-grown  crystals  is  shown  in 
Fig.  3(b).  It  corresponds  to  "Center  C"  in  Roelofs'  paper*  and,  thus,  in  our  notation  is 
Center  C^.  This  center  is  considerably  less  intense  than  the  previously  described  center. 
The  ENDOR  technique  shows  hyperfine  interactions  with  three  phosphorus  nuclei.  There 
were  no  ENDOR  lines  which  could  be  assigned  to  hydrogen. 

Figure  5  shows  the  infrared  absorption  spectra,  and  thus  the  relative  amounts  of 
OH“,  for  our  particular  crystals.  There  is  approximately  an  order  of  magnitude  more 
hydrogen  in  the  hydrothermally  grown  crystal. 
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C.  Conclusions 

We  have  shown  that  a  variety  of  Ti^  centers  can  be  created  at  77  K  in  KTP. 
Furthermore,  the  specific  centers  formed  in  a  particular  crystal  depend  on  which  tech¬ 
nique  was  used  to  grow  the  crystal  (other  unidentified  factors,  of  course,  may  also  play  a 
role).  Hydrogen  is  adjacent  to  the  dominant  Ti^'^  centers  in  hydrothermally  grown 
crystals,  but  not  in  flux-grown  crystals.  This  is  consistent  with  the  reports  of  previous 
investigators  that  hydrothermally  grown  crystals  contain  significantly  more  hydrogen  than 
flux-grown  crystals.^*^  In  general,  we  expect  an  OH"  ion  will  provide  the  stabilizing 
influence  needed  to  trap  an  electron  on  an  adjacent  titanium  ion.  Why  this  does  not  occur 
in  the  flux-grown  crystals  is  a  major  question.  Also,  it  is  not  known  what  adjacent  entity 
stabilizes  the  electron  on  a  titanium  ion  in  the  flux-grown  crystals.  The  most  likely 
possibility  at  this  time  is  an  oxygen  vacancy.  Clearly,  more  detailed  ENDOR  studies  are 
needed  to  clarify  the  nature  and  role  of  the  stabilizing  entities  associated  with  each  Tp"^ 
center  in  KTP. 

The  observation  that  a  variety  of  Ti^  centers  may  exist  in  this  material  has  sig¬ 
nificant  ramifications  in  selecting  specific  KTP  crystals  for  high-power  laser  applications. 
These  Ti^"*"  centers,  and  their  associated  hole  traps,  will  each  have  a  characteristic 
thermal  stability.  For  example,  at  room  temperature,  many  of  these  defects  may  ther¬ 
mally  decay  in  milliseconds,  whereas  others  may  have  lifetimes  of  minutes  or  even  hours. 
Knowing  which  centers  will  form  during  exposure  to  an  intense  laser  beam  and  knowing 
the  thermal  stability  of  those  centers  may  allow  a  crystal's  susceptibility  to  gray  tracking 
to  be  predicted.  Beyond  simply  developing  selection  criteria  for  gray-track-resistant 
crystals,  it  is  reasonable  to  assume  that  future  improvements  in  the  power-handling 
capabilitity  of  KTP  crystals  will  focus  on  the  elimination,  or  neutralization,  of  those 
stabilizing  entities  that  allow  an  electron  to  be  trapped  on  a  titanium.  In  this  last  regard, 
it  is  important  to  realize  that  elimination,  or  neutralization,  of  hole  trapping  sites  may  be 
more  realistic,  and  just  as  effective,  as  focusing  on  the  elimination  of  electron  traps. 


m.  IDENTinCATlON  OF  A  RADIATION-INDUCED  HOLE  CENTER  IN  KTP 

Potassium  titanyl  phosphate  (KTi0P04),  better  known  as  KTP,  is  widely  used  to 
frequency-double  infrared  lasers  into  the  visible.  However,  point  defects,  with  their 
associated  optical  absorption  bands,  can  limit  the  power-handling  capability  of  these 
doubler  crystals.  Optical  damage  in  KTP  usually  manifests  itself  as  small  gray  tracks 
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produced  when  a  crystal  is  subjected  to  high-power,  high-repetition-rate  laser 
pulses.  The  defects  responsible  for  the  gray  tracking  have  not  been  identified,  nor 

has  the  mechanism  by  which  these  defects  are  formed.  Before  significant  progress  can  be 
made  in  improving  the  performance  of  KTP  doubler  crystals,  in-depth  studies  are  needed 
to  identify  and  characterize  the  various  point  defects  that  can  exist  in  the  material.  A 
better  understanding  of  the  mechanisms  in  KTP  by  which  point  defects  form,  stabilize, 
interact,  and  decay  will  then  lead  to  a  correlation  of  gray  tracking  with  specific  defects. 

Very  few  investigations  of  point  defects  in  KTP  have  been  reported  thus  far. 
Scripsick  et  al.^  observed  three  paramagnetic  radiation-induced  point  defects  in  hydro- 
thermally  grown  KTP.  These  defects  were  only  stable  below  room  temperature  and 
included  a  trapped-hole  center,  a  Ti^"^  center,  and  a  Pt^'^  center.  Roelofs,*  in  a  study  of 
flux-grown  material,  reported  four  paramagnetic  Ti^"*"  centers  produced  by  electric  fields 
or  by  annealing  in  a  hydrogen  atmosphere.  After  x-ray  irradiating  flux-grown  crystals  at 
room  temperature,  Andreev  and  Efimov^  found  an  aluminum-related  paramagnetic  hole 
center  and  the  four  Ti^'*'  centers  previously  reported  by  Roelofs.  Morris  and  cowork- 
ers^^'i^  used  ionic  conductivity  results  from  KTP  to  obtain  evidence  for  potassium  and 
oxygen  vacancies.  Other  investigations*^’**’^  have  focused  on  the  large  number  of  OH‘- 
related  infrared  absorption  bands  in  KTP  in  the  3400  to  3700  cm**  spectral  range.  In 
addition,  there  have  been  a  series  of  magnetic  resonance  studies  of  transition-metal-ion 
impurities  (i.e.,  Fe^'*',  Cr^'*',  and 

In  this  section,  we  describe  the  results  of  a  detailed  electron  paramagnetic 
resonance  (EPR)  and  electron-nuclear  double  resonance  (ENDOR)  investigation  of  the 
dominant  low- temperature,  radiation-induced  trappied-hole  center  in  KTP.  From  the 
angular  dependence  of  the  EPR  and  ENDOR  data,  we  show  that  this  defect  has  the  hole 
localized  on  an  oxygen  ion  between  two  titaniums  and  that  a  nearest-neighbor  potassium 
vacancy  serves  as  the  stabilizing  entity.  This  particular  defect  is  the  principal  hole  trap  in 
both  hydrothermally  grown  and  flux-grown  KTP  crystals.  Its  identification  provides 
spectroscopic  evidence,  through  the  EPR  technique,  for  the  existence  of  potassium 
vacancies  in  KTP.  Also,  this  defect  is  expected  to  play  a  significant  role  in  the  gray¬ 
tracking  phenomenon  because  of  its  relatively  long  lifetime  (i.e. ,  milliseconds)  at  room 
temperature. 
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A.  Crystal  Structure  of  KTP 

The  crystal  structure  of  KTP  was  initially  determined  by  Tordjman  et  al.^^  and 
later  verified  by  Thomas  et  al..^  It  is  orthorhombic  with  space  group  Pna2j  and  the 
symmetry  elements  are  two  glide  planes  and  one  screw  axis.  A  unit  cell  of  this  material 
contains  64  ions  and  has  lattice  parameters  a  =  12.819  A,  b  =  6.399  A,  and  c  =  10.584 
A.  The  first  glide  plane  consists  of  a  reflection  through  a  plane  perpendicular  to  the  a 
axis  with  subsequent  translations  of  1/2  the  lattice  parameter  along  the  b  and  c  axes.  The 
second  glide  plane  is  a  reflection  through  a  plane  perpendic-ular  to  the  b  axis  with  a 
translation  of  1/2  the  lattice  parameter  along  the  a  axis.  The  screw  axis  is  a  180°  rotation 
around  the  c  axis  with  a  translation  of  1/2  the  lattice  parameter  along  the  c  axis. 

Titanium  ions  are  found  in  two  crystallographically  inequivalent  sites  within  the 
KTP  structure.  There  also  are  two  distinct  sites  for  the  potassium  ions  and  two  for  the 
phosphorus  ions.  The  oxygen  ions  occupy  ten  crystallographically  inequivalent  sites. 
Together,  these  ions  represent  16  distinct  sites  in  the  crystal.  Their  lattice  positions  are 
given  in  Table  I,  and  the  positions  of  the  remaining  atoms  in  the  unit  cell  can  be  gener¬ 
ated  by  applying  the  symmetry  elements  of  the  crystal.  Specifically,  a  second  set  of  16 
ion  positions  can  be  generated  from  the  x,y,z  values  in  Table  I  by  applying  a  transfor¬ 
mation  representing  the  first  glide  plane.  The  appropriate  equations  are 


x* 

=  1/2  -  X. 

(la) 

y’ 

=  y+1/2, 

(lb) 

and  z' 

=  z-hl/2. 

(Ic) 

In  this  and  the  subsequent  transformation,  all  of  the  x,y,z  values  are  in  units  of  the  a,b,c 
lattice  parameters,  respectively.  A  third  set  of  16  ion  positions  can  be  generated  from  the 
x,y,z  values  in  Table  I  by  applying  a  transformation  representing  the  second  glide  plane. 
An  appropriate  set  of  equations  in  this  case  are 


X' 

=  x+1/2. 

(2a) 

y’ 

=  3/2  -  y, 

(2b) 

and  z' 

=  z. 

(2c) 

The  fourth  set  of  16  ion  positions  can  be  generated  by  sequentially  applying  both 
of  the  glide  plane  transformations  given  above  and  is  equivalent  to  applying  a 
transformation  representing  the  screw  axis  alone. 
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Table  I.  Ion  positions  are  given  in  units  of  the  a,  b,  and  c  lattice  parameters 
for  one-quarter  of  a  unit  cell  in  KTP.  The  remaining  48  ions  in  the  unit  cell 
are  obtained  from  these  by  application  of  the  crystal's  symmetry  elements. 
fTaken  from  Ref.  25,  except  for  a  different  choice  of  Ti(l)  and  OT(l).] 


Atom 

x/a 

y/b 

z/c 

K(l) 

0.37807 

0.7806 

0.6880 

K(2) 

0.10526 

0.6990 

0.9332 

Ti(l) 

0.37290 

0.5001 

0.9996 

Ti(2) 

0.24658 

0.2695 

0.74836 

P(l) 

0.49808 

0.3363 

0.7397 

P(2) 

0.18079 

0.5020 

0.4872 

0(1) 

0.4859 

0.4867 

0.8497 

0(2) 

0.5103 

0.4657 

0.6170 

0(3) 

0.4004 

0.1986 

0.7208 

0(4) 

0.5934 

0.1930 

0.7589 

OT(l) 

0.2752 

0.4653 

0.8561 

OT(2) 

0.2232 

0.0413 

0.6097 

0(5) 

0.1126 

0.3106 

0.4585 

0(6) 

0.1113 

0.6918 

0.5117 

0(7) 

0.2525 

0.5402 

0.3718 

0(8) 

0.2528 

0.4619 

0.6008 

B.  EPR  Results 

The  crystal  used  in  this  study  was  grown  by  the  hydrothermal  technique  at  Airtron 
(a  Division  of  Litton)  in  Morris  Plains,  New  Jersey.  Seed  temperature  during  growth 
was  550°C  and  nutrient  temperature  was  570°C.  Pressure  within  the  autoclave  was  1.72 
X  10*  Pa  (i.e.,  25,000  psi)  and  growth  time  was  approximately  40  days.  EPR  and 
ENDOR  data  were  taken  on  a  sample  3x3x4  mm^  which  had  been  cut  from  the  larger 
as-grown  boule.  The  faces  of  this  rectangular  sample  were  perpendicular  to  the  three 
high  symmetry  directions  (i.e.,  the  a,  b,  and  c  axes). 
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The  KTP  sample  was  irradiated  at  77  K  with  x-rays  to  produce  the  trapped-hole 
centers.  Irradiations  lasted  up  to  four  hours  with  the  x-ray  tube  operating  at  60  kV  and 
30  mA.  Following  an  irradiation,  the  sample  was  transferred  without  warming  into  an 
Oxford  Instruments  Model  ESR-900  helium  gas  flow  system  for  the  subsequent  EPR  or 
ENDOR  measurements.  The  EPR  data  were  taken  with  a  Bruker  Model  ESP  300 
spectrometer  operating  at  9.440  GHz  with  100-kHz  static  field  modulation.  A  Vanan  E- 
500  digital  gaussmeter  was  used  to  measure  the  magnetic  field  and  a  Hewlett-Packard 
5340A  counter  was  used  to  measure  the  microwave  frequency.  A  small  MgO:Cr  crystal 
was  used  to  correct  for  the  difference  in  magnetic  field  between  the  KTP  sample  and  the 
gaussmeter  probe  (the  isotropic  g  value  for  Cr^"^  in  MgO  is  1.98(X)).  The  Bruker  ESP 
300  spectrometer  also  was  used  for  the  ENDOR  measurements.  In  these  experiments,  the 
rf  field  was  frequency  modulated  at  12.5  kHz  with  the  depth  of  modulation  varying  from 
10  to  50  kHz.  The  ENDOR  coil  was  helical  and  was  attached  to  the  Oxford  Instruments 
glassware  extending  through  the  Bruker  ENDOR  cavity. 

Before  irradiation,  the  KTP  sample  was  visibly  clear;  after  irradiation  at  77  K,  it 
had  a  reddish-brown  coloration.  Correspondingly,  only  a  weak  Fe^+  EPR  spectrum  was 
present  before  irradiation,  and  a  very  large  holelike  spectrum  was  present  after  the  77-K 
irradiation.  Other  EPR  spectra,  due  to  Ti^"^  and  Pt^"^  centers,  were  present  after  the 
irradiation  but  they  are  not  the  focus  of  this  investigation. 

The  trapped-hole  center's  EPR  spectrum,  taken  with  the  magnetic  field  parallel  to 
the  c  axis  of  the  crystal,  is  shown  in  Fig.  6(a).  It  consists  of  a  primary  doublet  having  a 
2.57-G  separation  and  secondary  triplets  (with  1:2:1  relative  intensities)  having  0.69-G 
separations.  The  g-value  along  the  c  direction  is  2.03(X).  Based  on  this  c-axis  spectrum, 
Scripsick  et  al.^  concluded  that  this  was  an  S  =  1/2  defect  interacting  with  three  I  =  1/2, 
100%  abundant  nuclei.  When  the  magnetic  field  is  along  the  a  direction,  the  spectrum 
has  a  g-value  of  2.0131  and  consists  of  a  doublet  split  by  about  1.5  G.  The  b-direction 
spectrum  has  a  g-value  of  2.0251  and  consists  of  a  doublet  having  a  splitting  of  2.54  G. 
These  latter  two  spectra  are  shown  in  Figs.  6(b)  and  6(c). 

Although  the  trapped-hole  center's  EPR  spectrum  is  best  observed  near  30  K,  it 
can  be  monitored  at  temperatures  above  80  K.  At  the  higher  temperatures,  the  triplet 
structure  observed  in  Fig.  6(a)  is  unresolved.  Also,  below  20  K,  there  is  a  broadening 
which  decreases  the  resolution.  If  the  temperature  of  the  sample  rises  above  approxi- 
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mately  160  K,  the  trapped-hole  center  anneals  and  the  sample  must  be  re-irradiated  before 
continuing  with  a  study  of  the  defect. 

For  an  arbitrary  direction  of  the  magnetic  field,  there  are  four  magnetically 
inequivalent  orientations  of  the  trapped  hole  center  and  there  are  four  sets  of  EPR  lines. 

If  the  magnetic  field  is  restricted  to  the  a-c,  b-c,  or  a-b  planes,  the  four  sites  are  pairwise 
degenerate  and  only  two  sets  of  EPR  lines  are  observed.  Finally,  when  the  magnetic 
field  is  along  any  of  the  three  high-symmetry  axes,  all  four  sites  are  magnetically  equiv¬ 
alent  and  only  a  single  set  of  EPR  lines  is  observed  (as  shown  in  Fig.  6).  In  the  present 
investigation,  EPR  data  were  taken  at  approximately  5°-intervals  in  the  a-c,  b-c,  and  a-b 
planes  of  the  crystal.  The  observed  angular  dependence  associated  with  the  g-matrix  is 
shown  in  Fig.  7  where  hyperfine  splittings  are  ignored  and  the  data  points  represent  the 
center  of  each  set  of  hyperfine-split  lines.  Although  one  hyperfine  interaction  was  re¬ 
solved  at  a  majority  of  the  angles  for  which  EPR  data  were  taken,  its  analysis  is  not 
attempted  until  the  ENDOR  data  is  available  later  in  this  section. 

The  spin-Hamiltonian  for  an  S  =  l/2, 1  =  1/2  defect  is 

H,  =  pS-g-H  -h  S  A-I-gNpNH-I  (3) 

where  the  terms  represents  the  electron  Zeeman  interaction,  the  hyperfine  interaction,  and 
the  nuclear  Zeeman  interaction,  respectively.  A  minimum  of  six  parameters  (three 
principal  values  and  three  Euler  angles  specifying  the  directions  of  the  principal  axes)  are 
needed  to  describe  the  g  matrix  and  another  six  are  needed  for  each  hyperfine  matrix.  In 
this  paper,  the  Euler  angles  are  converted  into  (0,(p)  pairs  of  angles,  thus  making  it  easier 
to  describe  the  direction  of  each  principal  axis.  The  polar  angle  0  is  measured  relative  to 
the  +c  direction  and  the  azimuthal  angle  cp  is  measured  relative  to  the  -f-a  direction  in  the 
c  plane  with  positive  rotation  being  from  a  to  b. 

The  angular  dependence  of  the  EPR  spectrum  was  used  to  obtain  the  g-matrix 
parameters.  In  particular,  a  least-squares-fitting  program  repeatedly  diagonalized  a  2  x  2 
Hamiltonian  matrix  (representing  only  the  electron  Zeeman  term)  while  using  95  inde¬ 
pendent  measurements  of  EPR  line  positions  and  their  associated  microwave  frequencies 
as  input  data.  These  95  magnetic  fields  represented  the  mid-points  of  sets  of  hyperfine- 
split  lines  (i.e.,  the  data  points  in  Fig.  7).  Table  II  contains  the  resulting  "best-fit"  g- 
matrix  parameters.  The  average  deviation  between  the  measured  and  the  final  calculated 
line  positions  was  0.085  G. 
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Table  II.  Spin-Hamiltonian  parameters  for  the  trapped-hole  center  in  KTP, 
The  g  matrix  was  obtained  from  EPR  data  and  the  hyperfine  matrices  were 
obtained  from  ENDOR  data. 


Principal  directions 

Matrix 

Principal  values 

0 

g 

2.0008 

57.90 

156.80 

2.0250 

79.90 

6O.40 

2.0422 

34.10 

315.10 

Ai  (3lP) 

-6.809  MHz 

131.50 

91.60 

-8.125  MHz 

44.30 

66.70 

-3.848  MHz 

77.30 

0 

d 

A2  (31P) 

0.918  MHz 

72.30 

138.10 

0.547  MHz 

71.60 

234.20 

2.598  MHz 

26.00 

7.20 

A3  (3‘P) 

0.500  MHz 

136.90 

120.00 

-0.360  MHz 

132.60 

310.60 

2.280  MHz 

84.70 

35.70 

The  spin-Hamiltonian  parameters  given  in  Table  II  correspond  to  one  of  the  four 
crystallographically  equivalent  sites  occupied  by  the  trapped  hole  center  in  KTP .  All  four 
sites  have  the  same  principal  values,  but  the  change  in  directions  of  the  principal  axes 
when  going  from  one  site  to  another  must  reflect  the  symmetry  elements  of  the  lattice. 
The  sets  of  angles  needed  to  describe  a  principal-axis  direction  at  the  other  three  sites  are 
obtained  from  the  initial  set  as  follows: 


Site  1: 

(01, 9i) 

Site  2: 

(0,,  180o-(pi) 

Reflection  in  the  a  plane. 

Site  3: 

(0,,36Oo-(pi) 

Reflection  in  the  b  plane. 

Site  4: 

(0,.  1800  -K  (pi) 

Reflections  in  both  planes. 
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Note  that  a  right-hand  coordinate  system  does  not  survive  in  these  transformations  of  the 
dir^tions  of  principal  axes.  However,  either  sense  along  a  specified  line  describes 
equally  well  the  direction  associated  with  the  corresponding  principal  value  of  a  spin- 
Hamiltonian  matrix. 

C.  ENDOR  Results 

An  electron-nuclear  double  resonance  spectrum  from  the  trapped-hole  center  in 
KTP  is  shown  in  Fig.  8.  This  spectrum  was  taken  with  the  magnetic  field  along  the  c 
axis  and  it  contains  six  distinct  pairs  of  lines  centered  on  the  free  nuclear-resonance 
frequency  (v^)  of  phosphorus.  Each  of  the  observed  hyperfine  interactions  is  less  than 
the  free  nuclear-resonance  frequency;  thus,  the  ENDOR  spectrum  consist  of  pairs  of  lines 
separated  by  A  and  centered  on  The  two  ENDOR  lines  forming  the  most  widely  split 
pair  in  Fig,  8  are  separated  by  7.297  MHz  and  correspond  to  the  phosphorus  nucleus 
responsible  for  the  2.57-G  hyperfine  splitting  observed  in  the  c-axis  EPR  spectrum  (see 
Fig.  6(a)).  A  second  pair  of  phosphorus  ENDOR  lines  are  separated  by  2.232  MHz  and 
a  third  pair  of  phosphorus  lines  are  separated  by  1.503  MHz.  The  average  of  these  latter 
two  separations  matches  the  0.69-G  triplet  hyperfine  splitting  observed  in  Fig.  6(a). 
Except  for  an  isolated  line  at  3.188  MHz,  the  remaining  ENDOR  lines  in  Fig.  8  represent 
phosphorus  hyperfine  interactions  which  are  unresolved  in  the  EPR  spectra. 

ENDOR  spiectra  were  taken  with  the  magnetic  field  along  the  a,  b,  and  c 
directions,  and  also  at  selected  orientations  of  the  magnetic  field  within  each  of  the  a-c, 
b-c,  and  a-b  planes  of  the  crystal.  These  latter  measurements,  two  in  each  plane,  were 
made  at  angles  approximately  25°  from  each  of  the  two  high  symmetry  directions 
defining  the  plane.  This  gave  a  total  of  15  distinct  ENDOR  spectra,  one  each  for  the 
three  high  symmetry  directions  and  two  each  for  the  six  orientations  between  the  high 
symmetry  directions,  from  which  to  obtain  parameters  describing  the  various  phosphorus 
hyperfine  interactions. 

The  experimentally  measured  separations  of  pairs  of  ENDOR  lines  representing 
three  of  the  ^^P  hyperfine  interactions  are  plotted  in  Figs.  9  and  10.  These  three  sets  of 
ENDOR  data  were  separately  fit  using  the  spin-Hamiltonian  given  by  Eq.  (3).  Specifi¬ 
cally,  the  4  X  4  matrix  derived  from  the  spin-Hamiltonian  was  repeatedly  diagonalized 
while  the  six  parameters  describing  the  hyperfine  matrix  (three  principal  values  and  three 
Euler  angles)  were  varied.  This  procedure  continued  until  a  "best-fit"  set  of  hyperfine 
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parameters  were  obtained  for  each  of  three  nuclei.  These  results  are  listed  in  Table  II. 
The  average  deviation  between  the  measured  and  the  final  calculated  values  for  the 
separation  of  pairs  of  ENDOR  lines  was  7.5  kHz,  12.5  kHz,  and  10.0  kHz,  respectively, 
for  the  Aj,  A2,  and  A3  nuclei.  Three  additional  ENDOR  spectra  are  resolved  in  Fig. 
8,  but  they  were  not  analyzed  in  the  present  investigation. 

The  following  condition  was  used  to  choose  the  absolute  signs  for  the  principal 
values  of  the  Aj,  A2,  and  A3  hyperfine  matrices  in  Table  II.  When  a  hyperfine  matrix  is 
reduced  to  its  isotropic  and  anisotropic  parts  (as  is  done  in  the  next  subsection),  the 
anisotropic  parameter  representing  the  classic  dipole-dipole  interaction  must  have  the 
same  sign  as  the  nuclear  magnetic  moment.  This  condition  forces  all  three  principal 
values  of  the  A 1  matrix  to  be  negative  and  leads  to  the  interesting  result  wherein  the 
isotropic  parameter  for  Aj  has  an  "anomalous"  sign  opposite  to  that  of  the  nuclear 
magnetic  moment.  Similar  "anomalous"  signs  for  the  isotropic  hyperfine  parameter  have 
been  observed  for  trapped-hole  centers  in  other  oxides.^^  The  situation  is  simpler  for  the 
A2  and  A3  matrices  where  the  isotropic  and  anisotropic  parameters  have  the  same  sign. 

D.  Model  of  the  Trapped-Hole  Center 

The  g  matrix  reported  in  Table  II  has  a  small,  but  significant,  positive  shift  which 
strongly  suggests  that  the  center  is  holelike,  and  thus  must  be  a  defect  consisting  of  a 
missing  electron  from  an  oxygen.  A  hole  residing  on  an  oxygen  ion  (i.e.,  an  O'  ion)  in 
an  oxide  material  typically  gives  positive  g  shifts  in  the  range  of  0.02  to  0.08.  This  is 
easily  seen  by  considering  the  g  values  predicted  for  an  unpaired  electron  occupying  an 
oxygen  p  orbital  in  a  tetragonal  crystal  field. To  first  order,  these  are 

g„  =  ge 

and  =  ge  -  2X/A.  (4b) 

The  parameter  A  is  the  splitting  between  the  singlet  and  doublet  energy  levels  associated 
with  the  three  p  orbitals  and  X  is  the  spin-orbit  constant.  For  an  oxygen  ion,  X  is  -135 
cm'^  This  gives  values  of  0.5  to  1.5  eV  for  A  when  the  g  shifts  are  in  the  range  of  0.02 
to  0.08.  In  this  simple  model,  A  is  interpreted  to  be  the  energy  associated  with  the 
optical  transition  from  the  ground  state  (p^)  to  the  first  excited  state  (p^.Py)  of  Oie  O'  ion. 
However,  the  optical  absorption  peak  positions  observed  for  hole  centers  in  oxides  are 
usually  a  factor  of  two  larger  than  the  value  of  A  predicted  from  their  measured  g  shifts. 
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In  the  case  of  the  trapped-hole  center  in  KTP,  the  p  orbital  containing  the  unpaired 
electron  is  in  an  orthorhombic  crystal  field.  This  causes  all  three  principal  values  of  the  g 
matrix  to  be  unique,^  in  agreement  with  the  results  in  Table  II. 

The  hyperfine  matrices  describing  the  interactions  between  phosphorus  nuclei  and 
the  unpaired  electron  associated  with  the  hole  center  are  given  in  Table  II.  It  is  instruc¬ 
tive  to  decompose  each  of  these  hyperfine  matrices  into  an  isotropic  part  and  an  aniso¬ 
tropic  part  (i.e.,  A  =  al  -I-  B).  The  anisotropic  portion  B  can  further  be  represented  in 
its  princip>al-axis  system  by  two  interaction  constants  as  follows.^ 


In  this  form,  the  parameter  b’  represents  the  deviation  from  axial  symmetry  for  the 
anisotropic  interaction.  Values  for  these  three  hyperfine  constants  (deduced  from  the  A 
matrices  in  Table  II)  are  aj  =  -6.261  MHz,  bj  =  1.206  MHz,  and  b^’  =  0.658  MHz  for 
matrix  Ap  they  are  a2  =  1.354  MHz,  62  =  0.622  MHz,  and  62’  =0.185  MHz  for 
matrix  A2,  and  they  are  aj  =  0.807  MHz,  bj  =  0.737  MHz,  and  b3'  =  0.430  MHz  for 
matrix  A3. 

As  noted  earlier,  the  g  matrix  provides  evidence  that  the  trapped  hole  resides  on 
an  oxygen  ion.  It  remains  to  be  determined  which  of  the  ten  inequivalent  oxygen  ions  in 
the  unit  cell  will  trap  the  hole.  Fortunately,  the  isotropic  constant  (aj  =  -6.261  MHz) 
extracted  from  the  largest  phosphorus  hyperfine  matrix  allows  one  to  significantly  narrow 
the  number  of  possibilities.  This  isotropic  constant  is  an  order  of  magnitude  smaller  than 
the  value  expected  if  the  oxygen  ion  containing  the  hole  were  a  component  of  a  PO4 
group.  Supporting  evidence  for  this  statement  is  provided  by  gamma-irradiated  KH2PO4 
(KDP)  where  a  trapped-hole  center,  consisting  of  an  unpaired  electron  localized  on  an 
oxygen  ion  in  a  PO4  group,  has  an  observed  isotropic  phosphorus  hyperfine  interaction  of 
96.7  MHz.^®  Thus,  because  of  the  small  value  of  the  phosphorus  isotropic  interaction  in 
KTP,  the  eight  oxygen  ions  directly  bonded  to  phosphorus  ions  are  eliminated  as  sites  for 
the  trapped  hole  and  only  the  two  oxygen  ions  without  an  adjacent  phosphorus  neighbor 
are  left  as  candidates.  These  latter  two  oxygen  ions  have  two  titaniums  as  nearest 
neighbors. 
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Additional  evidence  that  places  the  trapped  hole  on  one  of  the  two  oxygen  ions 
that  are  bonded  only  to  titanium  nearest-neighbors  comes  from  the  anisotropic  part  of  the 
largest  hyperfine  interaction.  The  parameter  b  in  Eq.  (5)  can  be  approximated  reasonably 
well,  in  our  case,  by  a  simple  point  dipole-dipole  interaction.  In  this  approximation,  b  is 
related  to  the  dipole  separation  R  by  the  following  expression.^ 

b  =  (Po/87t)  g^p^ggPg  (2/r3)  (6) 

A  brief  calculation  shows  that  a  bj  value  of  1.206  MHz  corresponds  to  a  separation 
distance  of  2.98  A  between  the  unpaired  electron  and  the  phosphorus  nucleus.  This 
distance  is  quite  different  from  the  nominal  oxygen-phosphorus  separation  of  1.54  A 
found  in  the  PO4  groups  within  KTP.  In  contrast,  the  calculated  value  of  2.98  A  com¬ 
pares  favorably  with  the  actual  distance  (i.e.,  3.22  A  to  3.32  A)  from  titanium-bridging 
oxygen  ions  to  their  nearest-neighbor  phosphorus  ions. 

A  remaining  step  in  developing  a  model  for  the  trapped-hole  center  is  to  explain 
how  the  hole  is  stabilized  at  a  particular  oxygen  ion.  There  is  no  evidence  in  the  litera¬ 
ture  to  support  the  stable  self-trapping  of  a  hole  on  an  oxygen  ion  in  an  otherwise  perfect 
oxide  lattice.  The  self-trapping  of  holes  does  occur  in  the  alkali  halides,^^  but  not  in 
oxides.  Thus,  a  trapped-hole  center  in  an  oxide  crystal  includes  a  component  that 
provides  the  stabilizing  force  to  maintain  the  hole  at  the  specific  defect  site  and  prevent  it 
from  migrating  (i.e.,  hopping)  through  the  lattice.  In  many  oxides,  the  stabilizing  entity 
is  a  nearest  neighbor  cation  vacancy and  we  suggest  a  similar  explanation  for  the 
trapped-hole  center  in  KTP.  Specifically,  we  propose  that  a  potassium  vacancy  is 
adjacent  to  the  oxygen  ion  trapping  the  hole  in  KTP.  In  support  of  this  expectation,  we 
find  reasonable  agreement  between  the  direction  of  the  unique  axis  of  the  g  matrix  (i.e., 
the  2.0008  value)  and  the  direction  from  a  titanium-bridging  oxygen  to  a  potassium  (see 
Table  III). 

A  trapped-hole  center  consisting  of  an  O'  ion  between  two  titanium  ions  with  a 
nearest-neighbor  potassium  vacancy  could  exist  in  four  distinct  configurations  within  the 
KTP  lattice.  The  four  possibilities  arise  because  there  are  two  crystallographically 
inequivalent  oxygen  sites  between  titanium  ions  and  there  are  two  crystallographically 
inequivalent  potassium  sites.  However,  only  one  of  the  four  possibilities  provides  a  good 
match  to  the  directions  of  the  unique  axes  of  the  experimentally  measured  spin-Hamilton- 
ian  matrices.  This  is  shown  in  Table  III  where  the  directions  and  separations  of  pairs  of 
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ions  in  the  "perfect"  KTP  lattice  are  compared  with  the  unique  axes  of  the  g  matrix  and 
the  hyperfine  matrices. 

Significant  lattice  relaxation  in  and  around  the  defect  (on  the  order  of  3  to  6%)  is 
expected  to  accompany  the  formation  of  the  potassium  vacancy  and  O"  ion,  and  this  will 
cause  the  directions  of  the  unique  principal  axes  of  the  g  and  hyperfine  matrices  to 
deviate  by  small,  but  measurable,  amounts  from  the  ion-to-ion  directions  in  the  perfect 
lattice.  Also,  the  spatial  distribution  of  the  p-orbital  wave  function  occupied  by  the  hole 
will  cause  deviations  of  hyperfine  principal  axes  from  the  ion-to-ion  directions. 


Table  III.  Comparison  of  directions  and  ion  separations  in  the  KTP  crystal 
structure  with  experimental  directions  and  separations  obtained  from  the  g 
matrix  and  the  hyperfine  matrices. 


Ion  Pair^  or 

Principal  Axis 

0 

Direction 

<P 

Separation 

Distance^ 

OT(l)— -K(2) 

72.84° 

145.53° 

2.765  A 

g-matrix  unique  axis 

57.9° 

156.8° 

t 

1 

1 

1 

1 

t 

O 

H 

67.50° 

163.89° 

3.219  A 

Aj-matrix  unique  axis 

77.3° 

p 

O 

2.98  A 

P(2) - OT(l) 

17.52° 

349.02° 

4.094  A 

A2-matrix  unique  axis 

26.0° 

7.2° 

3.72  A 

P*(l) - OT(l) 

73.05° 

28.51° 

4.226  A 

A3-matrix  unique  axis 

84.7° 

35.7° 

3.52  A 

^  The  position  of  ion  P'*(l)  is  obtained  from  ion  P(l)  by  applying  the  second  glide  plane 
(i.e.,  Eqs.  2)  and  then  translating  one  unit  in  the  -a  direction  and  one  unit  in  the 
-b  direction. 

2  For  the  hyperfine  principal  axes,  this  separation  distance  is  the  classic  dipole-dipole 
prediction  as  given  by  Eq.  (6). 
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However,  despite  these  uncertainties  associated  with  lattice  relaxations  and  the  lack  of 
corrections  for  the  spatial  distribution  of  the  unpaired  electron's  wave  function,  the 
agreement  found  in  Table  III  provides  encouraging  support  for  the  following  model. 
Specifically,  we  propose  that  the  trapped-hole  center  in  KTP  consist  of  a  hole  (i.e.,  a 
missing  electron)  localized  on  the  OT(l)  oxygen  ion  with  a  potassium  vacancy  at  the 
nearest-neighbor  K(2)  site.  This  oxygen  ion  bridges  between  the  Ti(l)  and  Ti(2)  titanium 
ions,  with  Ti(l)  providing  the  "long"  bond  to  OT(l)  and  Ti(2)  providing  the  "short" 
bond.  The  largest  hyperfine  interaction  is  with  the  P(l)  phosphorus  ion,  which  agrees 
with  the  lobe  of  the  p-orbital  pointing  toward  this  ion.  A  schematic  representation  of  this 
model  for  the  trapped-hole  center  is  shown  in  Fig.  11.  The  A2  and  A3  hyperfine  inter¬ 
actions  are  with  the  P(2)  and  P*(l)  phosphorus  ions  (see  Table  III  for  a  description  of  P*) 
and  both  of  these  ions  are  joined  via  oxygen  ions  to  Ti(2).  Since  Ti(2)  represents  the 
"short"  bond  to  OT(l),  it  is  reasonable  to  expect  that  more  of  the  spin  density  will  be 
transferred  to  the  Ti(2)  side  of  the  defect  rather  than  to  the  "long"  bond  Ti(l)  side. 

The  specific  model  illustrated  in  Fig.  1 1  suggests  that  only  potassium  vacancies  on 
the  K(2)  site,  or  one  of  the  other  three  equivalent  sites  in  the  unit  cell,  will  trap  a  hole  on 
an  adjacent  oxygen  ion.  This  can  be  explained  either  by  (a)  assuming  that  potassium 
vacancies  in  KTP  occupy  only  the  K(2)  site  and  not  the  K(l)  site,  or  (b)  assuming  that 
potassium  vacancies  occupy  both  sites  but  holes  are  only  trappied  by  the  vacancies  on  the 
K(2)  sites.  Choice  (a)  has  significant  ramifications  for  the  defect  structure  of  KTP,  if  it  is 
shown  to  be  true  by  later  experiments.  Choice  (b)  is  equally  likely  to  be  operative  since 
the  separation  between  the  oxygen  ion  and  the  potassium  vacancy  is  an  important  param¬ 
eter  in  deciding  whether  the  resulting  hole  center  will  be  thermally  stable.  As  described 
earlier,  there  are  four  combinations  of  titanium-bridging  oxygens  and  potassium  sites  in 
the  KTP  lattice  and  their  separations  are  2.722  A,  2.765  A,  2.996  A,  and  3.057  A  in  the 
perfect  lattice.  The  model  in  Fig.  1 1  corresponds  to  the  2.765  A  separation.  Subsequent 
investigations  in  which  the  sample  is  irradiated  and  maintained  at  temperatures  near  that 
of  liquid  helium  may  reveal  additional  stable  configurations  of  trappied  hole  centers  in 
KTP. 


In  concluding  this  section,  we  note  that  a  paramagnetic  trapped-hole  center  has. 
been  produced  in  KTP  during  x-ray  irradiation  at  77  K.  The  defect  is  stable  up  to  160  K, 
but  its  EPR  spectrum  can  best  be  observed  near  30  K.  There  is  a  strong  hyperfine 
interaction  with  one  ^*P  nucleus  and  weaker  interactions  with  additional  ^*P  nuclei. 
Parameters  for  the  g  matrix  and  three  hyj)erfine  matrices  were  obtained  from  the  angular 
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dependence  of  the  EPR  and  ENDOR  spectra,  respectively.  The  g-shifts,  the  isotropic 
hyperfine  interactions,  and  the  anisotropic  hyperfine  interactions  suggested  a  model  in 
which  the  hole  is  localized  on  a  bridging  oxygen  ion  between  two  titanium  ions.  The 
hole  is  stabilized  at  this  oxygen  by  a  nearest-neighbor  potassium  vacancy. 

rv.  EPR  OF  A  NATIVE  ACCEPTOR  IN  AS-GROWN  ZnGeP2 

Zinc  germanium  phosphide  (ZnGeP2)  is  rapidly  becoming  an  important  frequency 
conversion  material  for  the  infrared  spectral  region. This  material  is  transparent 
from  near  0.7  to  1 1  and  has  a  large  nonlinear  optical  coefficient,  thus  making  it  an 
excellent  candidate  for  second  harmonic  generation  (SHG)  and  optical  parametric 
oscillator  (OPO)  applications.  Furthermore,  phase  matching  can  be  achieved  at  room 
temperature  across  the  entire  range  of  transparency.  In  the  past,  the  development  of 
ZnGeP2-based  optical  devices  has  been  limited  (1)  by  a  lack  of  large,  properly  oriented 
single  crystals  and  (2)  by  an  acceptor-related  absorption  band  in  the  near  infrared 
extending  from  approximately  0.7  fim  to  2.5  ^m.  It  is  only  recently,  through  improved 
crystal-growth  procedures  and  post-growth  anneals,  that  major  progress  has  been  made  in 
overcoming  both  limitations.^^  These  efforts  have  culminated  in  several  demonstrations 
of  2.05-Mm-pumped  Type-I  ZnGeP2  OPOs.  Budni  et  al.^^  reported  a  26%  overall 
conversion  efficiency  and  a  37%  slope  efficiency  at  an  average  continuous  power  output 
of  0.6  W.  Schunemann  et  al.^"^  obtained  a  1.6-W  average  output  power  while 
maintaining  an  overall  efficiency  of  18%. 

Additional  improvements  in  the  performance  of  ZnGeP2-based  OPOs  can  be 
expected  if  the  material's  extrinsic  absorption  near  the  pump  wavelength  at  2.05  /im  is 
further  reduced.  A  deep  acceptor,  recently  designated^*  as  ALl  (i.e.,  acceptor  level 
one),  is  responsible  for  this  near-edge  absorption.  The  ALl  center  is  present  in  all 
ZnGeP2  crystals  at  concentrations  on  the  order  of  lO’^  to  10^®  cm'^;  however,  in  the 
past,  optical  measurements  combined  with  various  post-growth  treatments  have  not 
provided  sufficient  information  to  establish  a  microscopic  model  for  this  defect.  We  have 
noted  that  an  intense  electron-paramagnetic-resonance  (EPR)  signal^^  from  an  acceptor 
also  is  present  in  all  as-grown  ZnGeP2  crystals.  It  is  our  expectation  that  the  EPR  signal 
and  the  ALl  center  have  a  common  origin,  since  they  have  similar  concentrations,  both 
are  acceptors,  and  both  can  be  reduced  in  intensity  by  comparable  amounts  when  a  crystal 
is  thermally  annealed.  By  focusing  on  this  EPR  signal,  we  hope  to  better  understand  the 
ALl  center. 
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In  this  section,  we  describe  the  results  of  a  study  of  the  EPR  signal  which  is  pre¬ 
sent  in  as-grown  crystals  of  ZnGeP2.  The  initial  observation  of  this  signal  was  reported 
by  KieP^  in  1974.  We  have  made  a  more  careful  analysis  of  the  EPR  spectrum's  angular 
dependence  and  have  found  precise  principal  values  and  principal  axis  directions  for  the  g 
matrix  of  this  center.  Our  conclusions,  which  differ  considerably  from  those  of  Kiel, 
restrict  the  most  likely  model  for  this  acceptor  to  two  choices,  either  a  zinc  vacancy 
(V^jj)  or  a  zinc  ion  on  a  germanium  site  (Zn^^).  The  acceptor,  in  either  case,  is  in  a 
single-negative  charge  state  and  the  unpaired  spin  (i.e.,  the  hole)  is  localized  on  a  pair  of 
adjacent  phosphorus  atoms.  The  defect’s  primary  symmetry  axis  is  oriented  along  [011]- 
type  directions,  thus  giving  it  four  crystallographically  equivalent  orientations  in  the  lat¬ 
tice.  In  contrast,  KieP^  misinterpreted  the  angular  dependence  of  the  EPR  spectrum  and 
erroneously  concluded  that  the  pair  of  phosphorus  atoms  sharing  the  hole  had  their  pri¬ 
mary  symmetry  axis  along  [1 10]-type  directions. 

ZnGeP2  has  the  tetragonal  chalcopyrite  crystal  structure'*®  where  the  c  direction  is 
denoted  as  [001]  and  the  two  equivalent  a  directions  are  labeled  [100]  and  [010].  The 
lattice  parameters,  a  =  5.46  A  and  c  =  10.71  A,  give  a  value  for  c/a  of  1.96.  Thus,  the 
crystal  can  be  viewed  as  having  an  approximate  2%  distortion  from  the  simpler  zinc- 
blende  structure.  In  ZnGeP2,  each  21n  ion  and  each  Ge  ion  has  four  P  neighbors,  and 
each  P  ion  has  two  Zn  and  two  Ge  neighbors.  The  GeP4  units  are  nearly  perfect  tetra- 
hedra,  whereas  the  ZnP4  units  form  slightly  distorted  tetrahedra.'*®  In  particular,  the  two 
P-Zn-P  angles  in  the  ZnP4  unit  are  108.5®  and  11 1.4®,  as  compared  with  the  perfect 
tetrahedral  bond  angle  of  109.47®,  whereas  the  two  Ge-P-Ge  angles  in  the  GeP4  unit  are 
109.3®  and  109.6®. 

A.  EPR  Results 

The  single  crystals  used  in  this  study  were  grown  at  Lx>ckheed  Sanders  by  the 
horizontal  gradient  freeze  technique. They  are  nearly  stoichiometric  but  may  be  slight¬ 
ly  Ge-rich  because  of  the  greater  volatility  of  the  zinc  at  the  high  growth  temperatures. 
EPR  samples,  in  the  form  of  cubes  approximately  2.5  mm  on  a  side,  were  cut  from  larger 
boules.  They  were  oriented  such  that  [100],  [010],  and  [001]  directions  emerged  from 
the  three  orthogonal  faces.  A  Bruker  ESP-300  spectrometer  operating  at  9.45  GHz  was 
used  to  obtain  the  EPR  data.  Microwave  frequencies  and  magnetic  fields  were  measured 
with  a  Hewlett  Packard  5340A  counter  and  a  Varian  E-500  digital  gaussmeter,  respective- 
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ly.  A  small  MgO:Cr  crystal  was  used  to  correct  for  the  difference  in  magnetic  field 
between  the  sample  and  the  gaussmeter  probe  (the  isotropic  g  value  for  Cr^"''  in  MgO  is 
1.9800).  An  Oxford  Instruments  ESR-900  helium  gas  flow  system  maintained  the  sample 
temperature  near  20  K  during  the  measurements. 

The  EPR  spectra  taken  along  various  high  symmetry  directions  are  shown  in  Fig. 
12.  Each  of  these  spectra  consists  of  one  or  more  similar  triplets  (i.e.,  sets  of  three  lines 
equally  split  and  having  1:2:1  intensity  ratios)  which  preliminary  electron-nuclear  double 
resonance  (ENDOR)  experiments  have  shown  to  be  due  to  equal  hyperfine  interactions 
with  a  pair  of  phosphorus  nuclei.  A  complete  ENDOR  study  of  this  defect  in  ZnGeP2  is 
in  progress  and  will  be  reported  in  a  separate  paper.  The  EPR  data  in  Fig.  12  can  be 
described  by  the  S  =  1/2  spin-Hamiltonian 

%  =  pS-g-H  -I-  S-A-Ii  -I-  S-A-12 

where  the  first  term  represents  the  electron-Zeeman  interaction  and  the  last  two  terms 
represent  the  phosphorus  hyperfine  interactions.  The  focus  in  this  letter  is  on  the 
electron-Zeeman  term  and  the  information  provided  by  the  g  matrix. 

As  shown  in  Fig.  12,  the  EPR  spectrum  from  this  defect  consists  of  a  single  triplet 
when  the  magnetic  field  is  parallel  to  either  the  [001]  or  the  [110]  directions.  In  contrast, 
when  the  magnetic  field  is  along  the  [Oil]  direction,  there  are  three  distinct  triplets  with 
the  middle  one  being  twice  as  intense  (i.e.,  doubly  degenerate).  Finally,  with  the  mag¬ 
netic  field  in  the  [100]  direction,  two  poorly  resolved  triplets  are  present,  each  doubly 
degenerate.  From  this,  we  deduce  that  the  defect  has  four  crystallographically  equivalent 
orientations  in  the  lattice.  All  four  orientations  are  magnetically  equivalent  when  the 
magnetic  field  is  along  the  [001],  [110],  or  [1 10]  directions,  but  this  degeneracy  is 
partially  lifted  when  the  magnetic  field  is  rotated  in  either  the  { 100}  or  the  {110}  sets  of 
planes.  If  the  magnetic  field  were  along  an  arbitrary  direction,  the  orientational  degen¬ 
eracies  would  be  entirely  removed  and  the  EPR  spectrum  would  contain  four  distinct 
triplets,  although  they  might  be  poorly  resolved  because  of  the  relatively  small  g  shifts. 

A  g  matrix  having  principal  axes  along  the  [01 1],  [  1 00],  and  [0 1 1]  directions  fits 
the  data  presented  in  Fig.  12.  Given  that  the  [100]-  and  [010]-type  directions  are  equal 
and  the  [(X)l]  direction  is  unique,  there  are  four  possible  orientations  of  this  g  matrix  in 
the  chalcopyrite  structure  and  they  correspond  to  the  four  crystallographically  equivalent 
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orientations  observed  in  the  data.  The  principal  g  values  associated  with  these  principal- 
axis  directions  were  obtained  as  follows.  It  was  noted  that  the  lowest-field  triplet  in  the 
[Oil]  EPR  spectrum  had  the  magnetic  field  parallel  to  a  principal  axis  and  that  the  highest 
field  triplet  in  this  same  EPR  spectrum  had  the  magnetic  field  along  another  of  the 
principal  axes.  The  corrected  magnetic  fields  corresponding  to  the  centers  of  these  two 
triplets  are  3256.0  G  and  3373.4  G  (for  a  microwave  frequency  of  9450.70  MHz).  These 
data  give  2.074  and  2.002  for  the  two  principal  g  values.  The  remaining  g  value  is  2.021 
and  is  obtained  from  the  [110]  EPR  spectrum  where  the  fourfold  degenerate  triplet  is 
located  at  3328.1  G  (for  a  microwave  frequency  of  9450.00  MHz).  The  effective  g  value 
for  this  latter  line  depends  on  all  three  principal  g  values;  however,  two  of  the  principal 
values  were  already  known  at  this  point.  These  results  for  the  g  matrix  are  summarized 
in  Table  IV.  It  is  important  to  note  that  the  g-shifts  associated  with  this  EPR  center  are 
positive,  thus  indicating  that  the  defect  is  an  acceptor. 


Table  IV.  Principal  values  and  principal  axes  directions  of  the  g  matrix  for  the 
four  crystallographically  equivalent  defect  sites  labeled  1  through  4  in  Fig.  13. 


Principal  Values 

Site  1 

Principal  Axes  Directions 

Site  2  Site  3  Site  4 

2.002 

[Oil] 

[101] 

[TOl] 

[Oil] 

2.021 

[TOO] 

[010] 

[OTO] 

[100] 

2.074 

[Oil] 

[TOl] 

[101] 

[011] 

B.  Model  of  the  Acceptor 

The  observed  hyperfine  interactions  (i.e.,  the  triplets)  require  that  the  unpaired 
spin  be  equally  shared  between  two  phosphorus  atoms,  presumably  nearest  neighbors. 
This  is  consistent  with  the  observation  of  four  crystallographically  equivalent  orientations 
of  the  g  matrix.  Specifically,  four  phosphorus  atoms  surround  each  cation  site  (i.e.,  zinc 
or  germanium)  in  the  chalcopyrite  structure,  and  four  pairs  of  these  phos-phorus  atoms 
having  the  same  symmetry  as  the  g  matrix  can  easily  be  identified.  These  four  pairs. 
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labeled  1  through  4,  are  illustrated  in  Fig.  13.  Having  the  defect  occupy  only  these  four 
crystallographically  equivalent  sites  is  not  surprising  since  a  close  examination  of  the 
chalcopyrite  structure  reveals  that  the  six  possible  pairs  of  phosphorus  atoms  within  the 
tetrahedron  shown  in  Fig.  13  naturally  separate  into  two  sets.^®  Four  of  the  pairs  have 
one  separation  distance  and  the  remaining  two  pairs  have  a  different  separation  distance. 
For  example,  in  a  GeP^  unit,  the  fourfold-degenerate  pairs  have  separation  distances  of 
3.792  A  and  the  twofold-degenerate  pairs  have  separation  distances  of  3.798  A.  The 
difference  between  the  two  types  of  pairs  is  even  greater  in  a  ZnP^  unit  where  the  four¬ 
fold-degenerate  separation  is  3.856  A  and  the  twofold-degenerate  separation  is  3.924  A. 

Preliminary  ENDOR  data  provide  no  evidence  for  the  association  of  an  impurity 
with  the  defect,  and  this  is  consistent  with  the  growth  of  the  crystals  from  ultra-high- 
purity  (six  9s  pure)  starting  materials.  Furthermore,  impurities  are  not  expected  at  the 
10^^  to  10^®  cm*^  levels  exhibited  by  the  EPR  signal.  This  restricts  the  possible  models 
for  our  paramagnetic  center  to  native  defects  such  as  vacancies  and  antisite  atoms. 
Phosphorus  on  an  antisite  can  be  eliminated  from  consideration  because  the  observed 
hyperfine  is  too  small  by  more  than  an  order  of  magnitude.'*^  The  requirement  that  the 
unpaired  spin  be  localized  on  two  adjacent  phosphorus  atoms  eliminates  zinc  on  a 
phosphorus  site,  germanium  on  a  phosphorus  site,  and  a  phosphorus  vacancy. 

Germanium  on  a  zinc  site  is  eliminated  because  it  is  a  donor.  This  leaves  a  zinc  on  a 
germanium  site,  a  zinc  vacancy,  and  a  germanium  vacancy  as  the  possible  models  for  the 
negative-charged  acceptor.  The  germanium  vacancy  is  an  unlikely  candidate  because  the 
germanium-rich  nature  of  ZnGeP2  crystals  suggests  that  germanium  vacancies  are  not  a 
dominant  component  of  the  defect  structure  in  this  material.  Also,  the  germanium 
vacancy  is  an  unlikely  candidate  because  the  two  separation  distances  for  adjacent 
phosphorus  pairs  (where  the  unpaired  spin  resides)  are  only  slightly  different  which 
makes  it  likely  that  all  six  possible  pairs  of  phosphorus  atoms  surrounding  the  vacancy 
would  be  occupied,  in  disagreement  with  observation.  However,  lattice  relaxation  may 
invalidate  this  latter  argument.  In  conclusion,  we  are  left  with  either  a  zinc  vacancy  or  a 
zinc  on  a  germanium  site  as  the  model  for  our  EPR-active  acceptor  defect. 

ENDOR  experiments,  focused  on  the  germanium  and  zinc  neighbors,  are  the  only 
way  to  determine  which  of  the  two  proposed  models  is  correct.  Until  then,  both  models 
remain  plausible.  It  is  quite  possible  that  both  types  of  defects  may  be  present  in  the  as- 
grown  crystals  even  though  we  are  observing  only  one  in  a  paramagnetic  form.  The 
higher  volatility  of  zinc  during  crystal  growth  could  lead  to  large  concentrations  of  zinc 
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vacancies.  The  presence  of  large  numbers  of  zinc  and  germanium  antisites  also  is 
expected  since  the  material  solidifies  at  1020°C  as  a  zinc-blende  random  alloy  and 
undergoes  a  phase  transition  to  the  chalcopyrite  structure  at  952°C.  Zinc  and  germanium 
ions  are  randomly  distributed  over  the  cation  sites  in  the  high-temperature  phase,  and  not 
all  of  the  disorder  is  removed  when  the  crystal  is  cooled  below  the  phase  change.  Thus, 
large  concentrations  of  zinc  and  germanium  ions  may  exist  on  the  wrong  cation  site  even 
at  room  temperature.  Even  though  we  suggest  that  germanium  antisites  may  very  well 
exist  in  the  material,  we  believe  it  is  unlikely  they  are  directly  associated  with  the  EPR- 
active  defect  being  discussed  at  the  present  time. 

In  summary,  our  model  for  the  EPR-active  acceptor  defect  present  in  as-grown 
ZnGeP2  is  an  unpaired  spin  (i.e.,  a  hole)  shared  equally  by  a  pair  of  phosphorus  ions 
adjacent  to  either  a  zinc  vacancy  or  a  zinc  on  a  germanium  site.  The  EPR  experiments  do 
not  distinguish  between  the  two  possibilities.  Only  the  four  pairs  of  phosphorus  ions 
having  the  shorter  separation  and  oriented  along  [01  l]-type  directions  are  occupied  by  the 
hole.  This  model  is  illustrated  in  Fig.  13  where  the  four  crystallographically  equivalent 
orientations  of  the  acceptor  are  labeled  1  through  4.  At  this  time,  there  is  no  data  to 
suggest  whether  or  not  the  unpaired  spin  is  hopping  between  different  pairs  of  phosphorus 
atoms  within  a  given  tetrahedron.  Also,  there  is  no  evidence  concerning  the  identity  of 
other  defects  that  may  be  present  in  the  crystal  for  compensation  or  to  ensure  overall 
charge  neutrality. 

Thermally  annealing  an  as-grown  ZnGeP2  crystal  significantly  reduces'*^  the 
concentration  of  the  ALl  center  and,  similarly,  we  have  found  that  the  annealing  process 
also  decreases  the  concentration  of  the  EPR-active  acceptor.  Thus,  it  is  tempting  to 
suggest  that  the  ALl  center  and  the  acceptor  responsible  for  the  EPR  signal  are  the  same. 
Future  studies,  in  which  optical  absorption  and  EPR  measurements  are  made  on  the  same 
samples  after  irradiation  and  thermal  annealing  treatments,  wilt  aid  in  confirming  our 
tentative  identification. 

V.  EPR  SPECTRA  IN  AS-GROVVN  CdGeAs2 

Cadmium  germanium  arsenide  (CdGeAs2)  belongs  to  the  ternary  chalcopyrite 
family^®  of  nonlinear  optical  materials.  It  is  a  II-IV-V2  semiconducting  compound  analo¬ 
gous  to  the  widely  studied  III-V  materials  except  that  its  crystal  structure  is  tetragonally 
distorted  from  the  simpler  zinc-blende  structure  of  the  III-V  compounds.  At  room  tem- 
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perature,  the  lattice  parameters'*^  of  CdGeAs2  are  a  =  5.945  A  and  c  =  11.221  A  and 
the  space  group  for  the  crystal  is  M2d .  The  absorption  edge'*'*"'*^  of  CdGeAs2  occurs 
near  2.34  pm  (i.e.,  0.53  eV)  and  its  transparency  range  extends  to  17  pm. 

The  chief  interest  in  CdGeAs2  arises  from  its  extremely  large  nonlinear  optical 
coefficient.  The  reported  value'*^  of  35 1  pm/ V  for  d35  is  the  highest  among  known  com¬ 
pounds,  three  times  higher  than  that  of  ZnGeP2  and  nine  times  higher  than  AgGaSe2  (two 
commercially  available  materials  for  infrared  frequency  conversion  applications).  In 
addition,  CdGeAs2  has  sufficient  birefrigence  to  allow  phase-matching  throughout  its 
transparency  range.  Although  its  small  band  gap  precludes  the  use  of  CdGeAs2  for  most 
optical-parametric-oscillator  applications,  the  material  shows  promise  for  frequency¬ 
doubling  the  output  of  CO2  lasers.  For  example,  Iseler  et  al.'**  achieved  up  to  27%  con¬ 
version  efficiency  and  produced  almost  2  W  of  average  second-harmonic  output  power 
when  CdGeAs2  was  pumped  with  a  Q-switched,  10.6-pm  laser  operating  at  a  repetition 
rate  of  20  kHz.  Unfortunately,  these  investigators  had  to  cool  the  crystal  to  77  K  to 
reduce  the  optical  absorption  in  the  mid-infrared. 

A  major  challenge,  therefore,  in  the  development  of  CdGeAs2  crystals  for  second 
harmonic  generation  is  the  elimination  or  reduction  of  defect-related  absorption  bands  in 
the  region  of  the  output  wavelength.  Isomura  et  al.,'*'*  for  example,  describe  an  intense 
absorption  shoulder  extending  out  from  the  band  edge  in  a  p-type  crystal  having  a  hole 
concentration  of  2  x  10*^  cm'^.  In  a  more  detailed  study  characterizing  the  extrinsic 
absorption  in  the  2.5  to  17-pm  region,  Iseler  et  al.'*®  found  that  p-type  samples  could  be 
described  by  a  model  consisting  of  donors  (D)  of  undetermined  ionization  energy,  shal¬ 
low  acceptors  (Al)  that  are  closely  compensated  by  the  donors,  and  deep  acceptors  (A2) 
located  near  the  center  of  the  energy  gap.  They  observed  concentrations  of  the  A2  accep¬ 
tor  ranging  from  5  x  10*^  to  2  x  10^®  cm'^  in  samples  from  different  growth  runs,  thus 
leading  them  to  conclude  that  this  deep  acceptor  is  probably  a  native  defect,  rather  than 
an  impurity.^® 

The  dominant  acceptors  (and  donors)  in  as-grown  material  are  thought  to  come 
from  a  combination  of  (1)  nonstoichiometry  introduced  during  growth  and  (2)  the  incom¬ 
plete  removal  of  disorder  on  the  cation  sublattice  during  subsequent  cool-down.  This 
second  mode  of  defect  formation  is  expected  because  CdGeAs2  crystals  undergo  a 
disorder-to-order  phase  transition^*^  upon  cooling  through  630°C.  In  this  transition,  the 
crystal  goes  from  the  zinc-blende  structure  (where  Cd  and  Ge  randomly  fill  cation  sites) 
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to  the  chalcopyrite  structure  (where  Cd  and  Ge  occupy  alternating  cation  sites),  and  it  is 
quite  possible  that  some  residual  disorder  survives  to  room  temperature.  Thus,  the  opti¬ 
cally  active  point  defects  in  CdGeAs2  are  most  likely  various  types  of  vacancies  and 
antisite  ions.  Spectroscopic  techniques  capable  of  probing  the  bulk  of  the  crystal,  such  as 
photoluminescence  and  electron  paramagnetic  resonance  (EPR),  can  be  used  to  identify 
and  monitor  these  major  point  defects  in  CdGeAs2. 

In  this  section,  we  describe  the  results  of  an  EPR  investigation  of  as-grown 
CdGeAs2.  The  sample  used  in  this  study  was  cut  from  a  large  (19-mm  diameter  by  100- 
mm  length)  substantially  crack-free  single  crystal  grown  at  Lockheed  Sanders  by  the 
horizontal  gradient  freeze  technique.  Dimensions  of  the  EPR  sample  were  2.5  x  2.9  x 
2.7  mm^  along  the  [100],  [010],  and  [001]  directions,  respectively.  In  this  tetragonal 
lattice,  the  c  axis  is  defined  to  be  the  [001]  direction  and  the  two  equivalent  a  axes  then 
correspond  to  the  [100]  and  [010]  directions.  A  Bruker  ESP-300  spectrometer  operating 
at  9.45  GHz  was  used  to  obtain  the  EPR  data.  Microwave  frequencies  were  measured 
with  a  Hewlett  Packard  5340A  counter  and  magnetic  fields  were  measured  with  a  Varian 
E-500  digital  gaussmeter.  A  small  MgO:Cr  crystal  was  used  to  correct  for  the  difference 
in  magnetic  field  between  the  sample  and  the  gaussmeter  probe  (the  isotropic  g  value  for 
Cr^'^  in  MgO  is  1.9800).  An  Oxford  Instruments  ESR-900  helium  gas  flow  system 
maintained  the  sample  temperature  at  selected  values  between  8  and  40  K  during  the 
measurements.  Special  care  was  taken  in  cooling  and  in  warming  the  crystal  because  of 
the  significant  anisotropy  in  its  thermal  expansion  properties.  The  coefficient  of  thermal 
expansion  is  reported  to  be  ten  limes  greater  for  the  a  axis  than  the  c  axis.'*^ 

There  is  a  broad  EPR  signal  in  as-grown  CdGeAs2  which  can  easily  be  observed 
when  the  sample  temperature  is  reduced  below  60  K.  At  higher  temperatures,  the  Q  of 
the  microwave  cavity  is  destroyed  as  a  result  of  absorption  associated  with  free  carriers 
present  in  the  conduction  or  valence  bands  of  the  crystal.  The  EPR  spectrum  shown  in 
Fig.  14(a)  was  taken  with  the  magnetic  field  along  the  c  axis  of  the  crystal  and  the  spec¬ 
trum  shown  in  Fig.  14(b)  was  taken  with  the  magnetic  field  along  the  a  axis.  Spectro¬ 
meter  settings  were  the  same  for  both  traces;  the  sample  temperature  was  17  K,  the 
microwave  power  incident  on  the  cavity  was  7  mW,  and  the  amplitude  of  the  100-kHz 
modulation  was  10  gauss. 

In  both  of  the  spectra  shown  in  Fig.  14,  an  intense  EPR  signal  representing 
approximately  5  x  10*^  spins/cm^  is  observed  near  g  =  2.  The  overall  linewidth  in  both 
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traces  is  about  165  gauss,  with  structure  being  present  in  Fig.  14(a)  but  not  in  Fig.  14(b). 
There  was  no  change  in  the  linewidth  of  the  c-axis  spectrum  as  the  temperature  was 
reduced  from  40  K  to  8  K;  thus,  we  attribute  the  observed  linewidth  primarily  to 
unresolved  hyperfine  interactions  with  the  magnetic  nuclei  surrounding  the  paramagnetic 
defects.  We  cannot,  however,  rule  out  contributions  to  the  linewidth  from  random  strain- 
induced  variations  either  in  the  g  values  or  in  the  zero-field  splittings  (if  the  spin  S  is 
greater  than  1/2),  but  we  consider  these  latter  possibilities  to  be  less  likely.  If  hyperfine 
is  responsible  for  the  observed  linewidths,  a  significant  contribution  would  be  expected 
from  ’’'^As  nuclei  (I  =  3/2)  as  a  result  of  their  100%  abundance  and  lesser  contributions 
would  be  expected  from  ^^*Cd  and  **^Cd  nuclei  (both  I  =  1/2)  with  12.8%  and  12.2% 
abundance,  respectively.  Thus  far,  we  have  been  unsuccessful  in  obtaining  an  electron- 
nuclear  double  resonance  (ENDOR)  spectrum  from  the  EPR  signals  shown  in  Fig.  14.  In 
general,  the  ENDOR  experiment  requires  that  an  EPR  transition  be  saturated  (i.e., 
populations  equalized)  with  microwave  power.  This,  in  turn,  allows  nuclear  spin  flips  to 
be  detected  via  their  effect  on  the  saturated  EPR  signal.  The  ENDOR  experiments  in 
CdGeAs2  were  limited  by  our  inability  to  saturate  the  EPR  signal,  even  at  5  K, 
presumably  because  of  short  spin-lattice  relaxation  times. 

Close  examination  of  the  structure  present  on  the  c-axis  spectrum  in  Fig.  14(a) 
shows  it  to  be  slightly  asymmetric.  This  asymmetry,  which  is  even  more  apparent  in 
spectra  taken  at  lower  microwave  powers  and  higher  temperatures,  suggests  that  there  are 
at  least  two  distinct  spectra  (i.e.,  two  different  centers)  present  in  Fig.  14(a).  In  the  sim¬ 
plest  interpretation,  the  c-axis  spectrum  is  described  as  a  sharper  line  approximately  56 
gauss  in  width  superimposed  on  a  broad  line  165  gauss  in  width.  The  sharper  line  is 
slightly  offset  to  higher  field  and  has  a  g^.  value  of  2.015  while  the  broader  line  has  a  g^. 
value  of  2.018.  When  the  magnetic  field  is  along  the  a  axis,  as  in  Fig.  14(b),  the  spec¬ 
trum  consist  of  a  single  structureless  broad  line  with  g^  =  1.995.  Presumably,  this  a-axis 
spectrum  contains  unresolved  components  associated  with  each  of  the  signals  seen  in  the 
c-axis  spectrum.  During  rotation  of  the  magnetic  field  from  the  c-axis  to  the  a-axis,  i.e., 
going  from  Fig.  14(a)  to  Fig.  14(b),  the  shift  of  the  spectrum  occurs  gradually  and  no  in¬ 
creased  resolution  appears  at  intermediate  angles. 

Figure  15  shows  the  c-axis  EPR  spectrum  taken  at  41  K  with  a  reduced 
microwave  power  of  2.3  mW.  When  compared  to  Fig.  14(a),  this  higher-temperature 
spectrum  has  a  considerably  better  resolved  structure  but  a  much  worse  signal-to-noise 
ratio.  The  broad  underlying  line  which  was  dominant  in  Fig.  14(a)  is  still  observable 
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near  g  =  2,  but  now  there  are  additional  sharper  lines  superimposed  on  it  and  there  are 
small  outer  satellite  lines  near  2950,  3100,  and  3600  gauss.  Weak  intensities  and 
overlapping  components  prevented  us  from  following  the  angular  dependences  of  these 
sharper  central  and  satellite  lines.  It  is  not  yet  clear  whether  these  sharper  lines  arise 
from  a  single  additional  defect  or  they  represent  several  additional  paramagnetic  defects. 
Also,  a  preliminary  examina-tion  of  intensity  ratios  and  splittings  does  not  support 
assigning  the  satellite  lines  to  **^Cd  and  ^^^Cd  hyperfme  interactions. 

Our  initial  survey  of  EPR  spectra  in  as-grown  CdGeAs2  has  provided  evidence  for 
the  presence  of  several  paramagnetic  point  defects.  However,  there  is  not  sufficient 
resolution  in  the  spectra  to  conclusively  identify  the  specific  models  of  the  responsible 
defects.  Based  on  recent  work  in  other  ternary  chalcopyrite  crystals,  we  can,  neverthe¬ 
less,  tentatively  suggest  that  the  primary  EPR-active  defects  in  CdGeAs2  are  associated 
with  the  cation  sublattice.  Recently,  Rakowsky  et  al.^®  described  an  EPR  spectrum  from 
an  acceptor  in  as-grown  ZnGeP2  and  suggested  that  it  is  either  a  zinc  vacancy  (V2n)  or  a 
zinc  ion  on  a  germanium  site  (Zn^^).  The  unpaired  spin  in  ZnGeP2  interacts  primarily 
with  only  two  phosphorus  nuclei  neighboring  the  cation  site  and  has  g  values  of  2.002, 
2.021,  and  2.074.  This  anisotropy  in  g  values  corresponds  to  a  Ag  of  0.072.  In  contrast, 
the  present  investigation  gives  a  Ag  of  only  0.023  for  the  dominant  defect(s)  in  CdGeAs2. 
The  smaller  anisotropy  in  g  values  observed  in  CdGeAs2  suggests  that  the  unpaired  spin 
may  be  shared  by  the  four  As  ions  surrounding  the  cation  site,  thus  averaging  out  the 
anisotropy  that  would  occur  if  the  spin  were  localized  on  only  one  or  two  of  the  As 
neighbors.  Also,  a  rapid  hopping,  or  tunneling,  of  the  spin  among  the  As  neighbors 
would  agree  with  the  fast  spin-lattice  relaxation  encountered  in  our  attempts  to  obtain 
ENDOR  spectra.  An  additional  example  of  a  delocalized  spin  surrounding  a  cation 
vacancy  in  a  ternary  chalcopyrite  is  provided  by  the  recent  work  of  Aksenov  et  al.^^ 

They  observed  a  photo-induced,  isotropic  EPR  signal  at  g  =  2.019  in  CUAIS2  and 
assigned  it  to  a  copper  vacancy  (Vqj)  with  the  hole  equally  shared  by  the  four  nearest- 
neighbor  sulfur  ions.  Bardeleben  et  al.^^  previously  had  observed  a  similar  copper- 
vacancy  defect  in  CuGaS2. 
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FIGURE  CAPTIONS 


Figure  1.  EPR  spectra  from  centers  in  x-ray-irradiated  hydrothermally  grown  KTP 
crystals.  Center  is  shown  in  trace  (a)  and  Center  is  shown  in  trace 
(b). 

Figure  2.  ENDOR  spectrum  obtained  from  Center  A;,y  with  the  magnetic  field  parallel 
to  the  b  direction. 

Figure  3.  EPR  spectra  from  Ti^  centers  in  x-ray-irradiated  flux-grown  KTP  crystals. 
Center  is  shown  in  trace  (a)  and  Center  C^  is  shown  in  trace  (b). 

Figure  4.  ENDOR  spectrum  obtained  from  Center  By^^  with  the  magnetic  field  parallel  to 
the  b  direction. 

Figure  5.  Infrared  absorption  spectra  comparing  the  OH”  concentrations  in 
hydrothermally  grown  and  flux-grown  KTP. 

Figure  6.  EPR  spectra  of  the  trapped-hole  center  in  KTP.  The  data  were  taken  at 

approximately  30  K  and  the  magnetic  field  was  aligned,  in  turn,  along  each  of 
the  crystal's  high  symmetry  directions. 

Figure  7.  Angular  variation  (g  matrix  portion  only)  of  the  EPR  spectrum  of  the  trapped- 
hole  center.  The  experimental  data  are  represented  by  the  circles  while  the 
solid  lines  were  calculated  using  the  g  matrix  in  Table  II. 

Figure  8.  ENDOR  spectrum  of  the  trapped-hole  center  in  KTP  taken  with  the  magnetic 
field  along  the  c  axis  of  the  crystal.  The  measurement  temperature  was 
approximately  30  K  and  the  microwave  frequency  was  9.6627  GHz.  Stick 
diagrams  illustrate  the  three  pairs  of  phosphorus  ENDOR  lines  corresponding 
to  the  resolved  splittings  in  the  c-axis  EPR  spectrum  shown  in  Fig.  1(a). 

Figure  9.  Angular  variation  associated  with  the  largest  ^^P  hyperfine  interaction.  The 
circles  represent  ENDOR  data  and  the  solid  lines  are  calculated  using  the  Aj 
hyperfine  matrix  in  Table  II. 
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Figure  10.  Angular  variation  associated  with  the  A2  (dashed  lines)  and  A3  (solid  lines) 
hyperfine  interactions.  The  circles  represent  experimental  data  points  and  the 
lines  are  calculated  using  the  "best-fit"  hyperfine  matrices  in  Table  II. 

Figure  11.  Model  of  the  trapped-hole  center  in  KTP.  A  projection  of  the  lattice  on  the 
a-c  plane  is  shown  with  the  height  of  each  ion  along  the  b  direction  given  in 
angstroms.  The  hole  is  localized  in  a  nonbonding  p  orbital  of  an  oxygen  ion 
bridging  two  titanium  ions.  The  p  orbital  points  toward  a  potassium  vacancy 
(represented  by  the  large  square)  on  one  side  and  toward  a  phosphorus  ion  on 
the  opposite  side. 

Figure  12.  EPR  spectra  taken  at  20  K  in  ZnGeP2  with  the  magnetic  field  parallel  to  the 
[001],  [110],  [100],  and  [011]  directions.  The  "stick  diagrams"  below  each 
trace  indicate  the  triplet  positions. 

Figure  13.  The  four  crystallographically  equivalent  orientations  of  the  EPR-active 

acceptor  in  ZnGeP2.  The  center  of  the  tetrahedron,  labeled  M,  is  suggested 
to  be  either  a  zinc  vacancy  or  a  zinc  on  a  germanium  site. 

Figure  14.  EPR  spectra  of  as-grown  CdGeAs2  taken  at  17  K  with  (a)  the  magnetic  field 
parallel  to  the  c  axis  and  (b)  the  magnetic  field  parallel  to  the  a  axis. 

Figure  15.  EPR  spectrum  of  as-grown  CdGeAs2  taken  at  41  K  with  the  magnetic  field 
parallel  to  the  c  axis. 
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